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We study the effect of the crowded nature of the cellular cytoplasm on the translocation of a 
polymer through a pore in a membrane. By systematically treating the entropic penalty due to 
crowding, we show that the translocation dynamics are significantly altered, leading to novel scaling 
behaviors of the translocation time in terms of chain length. We also observe new and qualitatively 
different translocation regimes depending upon the extent of crowding, transmembrane chemical 
potential asymmetry, and polymer length. 
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Transport of a variety of biopolymers across a dividing 
membrane is a fundamentally important process in bio- 
logical systems [l| . Examples include transfer of proteins 
across cellular membranes or endoplasmic reticulum Q, 
gene swapping through bacterial pili 01 and RNA trans- 
port through nuclear pore complexes [4|]. Technological 
applications include gene delivery [5(] and DNA sequenc- 
ing Q . A considerable amount of theoretical work has 
focused on both the basic physics HI U| underly- 

ing the translocation process and on how details such as 
polymer-pore interactions 1 2L in trinsic polymer struc- 
ture [IH and confinement [la] affect the dynamics of 
the process. One important aspect that has received very 
little attention is the effect of crowding on the transloca- 
tion dynamics. Crowding due to macromolecular aggre- 
gates and other inclusions in the cellular cytoplasm can 
be as high as 50% by volume [l(| and is known to have 
considerable influence on reaction rates, protein folding 
rates and equilibria in vivo 13, [18( . A polymer threading 
its way through such a crowded environment is subject to 
a large entropic penalty which should dramatically affect 
the translocation dynamics. In this Letter we present 
the first systematic study of polymer translocation in 
terms of crowding. We show that the free energy penalty 
due to crowding has terms that scale sublincarly with 
the polymer length. The presence of non-linear terms in 
the free energy leads to qualitatively different transloca- 
tion dynamics including novel power law scalings of the 
translocation time with polymer length as well as situa- 
tions where the translocation time is nearly independent 
of the polymer length over several orders of magnitude. 

We consider a Gaussian polymer of length N (in units 
of the Kuhn length b) threading itself through a pore in 
a dividing membrane from cis (left) side to trans (right) 
side, as illustrated in Fig. [TJ The pore is assumed to 
be small enough that it allows only one monomer to 
pass through at a time, with an effective diffusion con- 
stant Dp. Crowding is modeled by randomly distributed 
spherical obstacles, sterically interacting with the poly- 
mer, of radius a and diffusion constant D a at a volume 
fraction <p c (<p t ) on the cis (trans) side. There could also 



be an excess chemical potential difference for monomers 
between cis and trans side, A/i. We now assume that we 
can treat the process quasi-statically with the polymer 
segments on both sides of the membrane being in equilib- 
rium at all times. The validity of the assumption depends 
on the relative magnitudes of the three time scales in the 
problem: the total translocation time r ~ (b 2 /Dp)f, the 
polymer relaxation time Tp ~ (rjb 3 /kpT)N 2 [l9j], and the 
timescale set by obstacle motion r Q ~ q L r 2 l' i /D a . Here f 
is a dimensionless function that characterizes the translo- 
cation process and r) is the medium's viscosity. Assuming 
equilibrium statistics for the polymer segments necessar- 
ily requires that Tp <C t Q. In the presence of obstacles, 
this assumption remains valid in two different regimes, 
To Ctj <Ct and Tp <C r <C r D . The first regime cor- 
responds to the situation where the obstacles diffuse fast 
enough that the "polymer segment + obstacles" system 
can be assumed to be at equilibrium (dynamic obstacles). 
In the second regime, the obstacles are essentially immo- 
bile on the translocation timescale and the polymer seg- 
ment achieves equilibrium statistics in this static obstacle 
environment. It should be noted that if the pore friction 
is not high enough (i.e. Tp ~ r) [9] or if Tp < t < r the 
quasi-static assumption breaks down leading to anoma- 
lous dynamics. In the regimes where the assumption re- 
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FIG. 1: Schematic illustration of the translocation process of 
a polymer in the presence of crowding 
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mains valid we have a well defined free energy barrier 
whose form is governed by the polymer statistics, pres- 
ence of the membrane, chemical potential gradient and 
the presence of crowding. Since the contributions to the 
free energy from factors other than crowding have been 
worked out before @, @| , we focus on the entropic penalty 
that arises from crowding. 

Polymer configurations in the presence of static ob- 
stacles correspond to Brownian walk trajectories with a 
diffusion constant b 2 /6 that have survived to a time t = n 
with the obstacles playing the role of traps. The fraction 
of allowed polymer configurations is therefore identical to 
the survival probability of such a Brownian walker, lead- 
ing to the free energy expression for the entropic penalty 
in units of ksT as 



F s cr (n) = -log S{n), 



(1) 
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(1\ ^/5 n 3/5 for n>iVx , ( 2 ) 

where N x ~ (a/6) 2 -3 / 2 represents the crossover poly- 
mer length from the Smoluchowski regime to the DV 
regime [25J . For the dynamic obstacle case, the "polymer 
+ obstacles" system comes to equilibrium. The presence 
of the obstacles gives rise to a depletion induced attrac- 
tive interaction between chain monomers that can lead 
to a collapsed polymer phase similar to that induced by 
poor solvent conditions [2(|. Simulations 27 1 and ana- 
lytical work [28|, [29| have shown that hard spheres can 
cause polymers to collapse if the sphere density is high 
enough. For the range of <f> (0.1 < <fi < 0.5) that we are 
interested in, the collapse occurs whenever a > 2b. For 
our purposes it is then safe to assume that our Gaussian 
polymer is in a collapsed "dense globule" state with a vol- 
ume Nb 3 . The free energy penalty is then given by the 
sum of the confinement entropic penalty (= tt 2 Rq/6R 2 
where R ~ N 1 / 3 b is the confining radius) and the work 
done to create a cavity of volume Nb 3 that is devoid of 
obstacles^ which is known exactly from scaled particle 
theory [30(. The resulting free energy penalty for the 
chain due to dynamic obstacles reads as 
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where S(t) is the survival probability of an appropri- 
ate Brownian walker at time t [HI, For short 
times, S(t) is given by the Smoluchowski solution S(t) ~ 
cxp {—nt) [illli^ ]. while for long times it is dominated by 
walkers trapped in large void regions giving rise to the 
stretched exponential Donsker-Varadhan (DV) solution 
S(t) ~ exp (-At 3 / 5 ) (Hill], where k and A are constants 
that depend on trap radius, trap density and geometry. 
Using Eq. fl} and the exact solutions for S(t) [ll|,[23|], we 
can explicitly compute the free energy penalty for a chain 
of length n due to static obstacles at volume fraction 0, F cr (<j),ri) 
yielding 
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FIG. 2: Translocation time r (in units of b 2 /Dp) of a poly- 
mer of length N releasing out of the crowded cis side with 
static obstacles (solid lines: 4> c — 0.1,0.3,0.5 in descending 
order), dynamic obstacles (dotted lines: for the same values 
of 4> c ), and no crowding (dashed-dotted line: cj) c = 0). Inset: 
r vs. obstacle volume concentration 4> c for static (solid lines: 



N — 10 , 10 , 10 in descending order) and dynamic obstacles 
(dotted lines: for the same values of N). Note <f>t = A/i = 
here and the value of b/a = 0.3 is used throughout the paper. 




The general expression for total free energy, taking into 
account the presence of the dividing wall and chemical 
potential difference, of a chain with n monomers on the 
trans side and N — n on the cis side is then given by 



F' 



,N-n) + F* r (<p t ,n) 



+ - ln[n(N - n)] + nAfi 



(4) 



for the static obstacle case and 



Ff ot (n) = F d cr (0 C , N - n) + F c d r (cf> t , n) + nA^ (5) 

for the dynamic obstacle case. It is to be noted that the 
logarithmic term in Eq. ((4]) , resulting from reduced chain 
configurations confined in a half-space, does not appear 
in Eq.© because Ff r already takes into account the en- 
tropic penalty associated with confining the chain to a 
dense globule. The translocation process can now be de- 
scribed by diffusion along the translocation coordinate 
n in the presence of a well-defined free energy barrier, 
F^\n). The dynamics of this process is governed by 
a Fokker-Planck equation which then allows one to com- 
pute the translocation time (mean first passage time for 
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the chain to diffuse across the pore [7J) as 

r d(s) = ^f 'dn f dn> e"^V). (6) 

Since we have nonlinear terms in the free energy we first 
consider the behavior of r for a general free energy func- 
tional with a power law scaling, e.g., F(n) ~ n a . Eq. (J6j) 
and saddle point approximations (in the large TV limit) 
to do the integrals yield 

r - N 2 -° for F(n) ~ (TV - n) a 

r - exp(TV a ) for - n a (7) 

for translocation out of and into a crowded half-space, 
respectively. Thus we anticipate new exponents charac- 
terizing the scaling properties of the mean passage time 
with the number of monomers in the crowded medium. 

We now consider, in detail, the impact of crowding 
on translocation dynamics in some physically interesting 
situations. The first example is that of a polymer es- 
caping from a crowded environment, i.e., where the cis 
side is crowded and the trans side is not {<pt = 0). In 
the static case, the dominant contribution to the free en- 
ergy has the form (j) 2 J b (N — n) 3 / 5 for a long chain and 
4> C (N — n) for a short chain (see Eq. ((2])). The transloca- 
tion times obtained by saddle point approximation scale 
as 4> c 2 l°N 7 / b and 4>~ l N, respectively. Exact numerical 
evaluations of Eq. ([6]) confirm the predicted power law 
scalings with both TV (Fig. [2]) and 4> c (inset in Fig. [2]). 
In contrast to the dynamic case where the leading or- 
der behavior simply corresponds to having an effective 
"osmotic pressure" from obstacles, the driving force for 
translocation in randomly distributed immobile obstacles 
is weakened by the existence of rare, large voids on the 
cis side that sufficiently long polymers can explore. This 
is the physical origin for the novel exponents describing 
the scaling of r with respect to both TV and <fi c for a long 
enough chain in the static case. 

Another situation of significant practical interest is 
translocation into a crowded environment driven by a 
chemical potential gradient. Here we take <f> c — and 
A/i = — jji that favors translocation into the crowded 
trans side. For the static obstacle case in the large TV 
limit, the linear chemical potential term always domi- 
nates for any non-zero value of (i. Thus r ~ TV in this 
limit, while, for fi = 0, r ~ exp(TV 3 / 5 ) (from Eq. (0). 
The situation is similar for the dynamic case but only if fi 
exceeds a critical threshold fi* that is sufficient to over- 
come the osmotic pressure term (/i* is here defined by 
the value of the prefactor of the n-linear osmotic pressure 
term in F cr ). For ji < fj,*, r ~ exp(TV) for a long chain 
when obstacles are mobile. Fig. 02a) and (c) clearly show 
these distinct limiting behaviors. The plots also reveal a 
striking phenomenon that seems to occur at intermedi- 
ate length scales. Depending on the parameter values, 
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FIG. 3: Translocation time r (in units of b 2 /Dp) vs. TV 
for (a) static and (c) dynamic obstacles on the trans side 
(4>t = 0.3, <f)c = 0) for different values of chemical potential 
gradient, showing the asymptotic power law scalings as well 
as the plateau regime at intermediate length scales. Note u* 
denotes the value of the prefactor of the n-linear term of F cr 
in Eq. (T2J) and ©, respectively. Different scaling behaviors 
(exponential, plateau, power law) of r in the (cj>t,n) phase 
space for (b) static and (d) dynamic obstacles at a fixed TV = 
10 5 . 



there appear to be regimes spanning several orders of 
magnitude in polymer length, where the translocation 
time is nearly independent of TV. The reason for this can 
be understood by considering the form of the relevant 
free energy profile as a function of the translocation co- 
ordinate n. The situation is particularly simple for the 
static obstacle case where there is a competition between 
a linear term (~ n) due to the chemical potential and a 
sublinear term (~ n 3 / 5 ) that comes from the crowding 
for a long chain. This gives rise to a free energy barrier 
whose height and position are nearly independent of the 
total length TV followed by a "downward slope" all the 
way to n = TV. The time taken to surmount this barrier, 
which is independent of TV, is rate-limiting and hence ef- 
fectively the total translocation time for polymer lengths 
shorter than the value of TV at which the time taken to 
traverse the downward slope becomes comparable. At 
this point the scaling crosses over to being linear in TV. 
The situation for the dynamic obstacle case is similar ex- 
cept that the presence of two sublinear terms "softens" 
the plateau because barrier height and position are no 
longer independent of TV. Figs. EHJb) and (d) show the 
different translocation time scaling regimes in the (</>t , fi) 
phase space at a fixed TV. As one might expect, large val- 
ues of /i imply the dominance of the chemical potential, 
leading to r ~ TV, while low values give exponential be- 
havior. The disappearance of the plateau regime for the 
dynamic obstacles at high concentrations is because the 
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FIG. 4: Translocation time r (in units of fo 2 /Dp) asafunction 
of chain length JV for (a) static(<^t = 0.05) and (c) dynamic 
obstacles(<ji> t = 0.2) for different values of <f> c - Different scal- 
ing behaviors (exponential, plateau, power law) of r in the 
(4>t,(t>c) phase space for (b) static and (d) dynamic obstacles 
at a fixed N = 3 x 10 4 . No chemical potential gradient is 
applied (A/i = 0). 



presence of two sublinear terms becomes more apparent 
at higher (j) t . 

As a final example we consider the polymer translo- 
cation when both sides are crowded and Afj, = 0. For 
the dynamic obstacles, qualitative picture is similar to 
that with trans side crowding and a chemical potential 
gradient (compare Fig. [3]Jc,d) with Fig. BJc,d)). It is to 
be noted however that here the presence of sublinear free 
energy contributions from both sides of the membrane 
smears out the plateau even more. In the static case 
however, the absence of a linear term in the long polymer 
limit gives rise to qualitatively different results. As seen 
in Fig. Sfa), even a minute amount of <pt leads to expo- 
nential translocation times at large polymer lengths de- 
spite substantial crowding on the cis side. This counter- 
intuitive behavior arises because the free energy profile in 
this situation always has a barrier whose height scales to 
leading order as A 3 / 5 , which in turn implies exponential 
barrier crossing times. 

In conclusion we have presented a systematic study of 
the impact of crowding on translocation dynamics. We 
have shown how crowding can lead to the emergence of 
novel exponents characterizing the scaling of transloca- 
tion time with polymer length and crowding volume con- 
centration. The existence of regimes where the transloca- 
tion time depends very weakly on polymer length, apart 
from being of theoretical significance, suggests the pos- 
sibility of designing filters for a tunable range of poly- 
mer lengths and also has implications for "synchronized" 
transport of proteins or nucleic acids from a wide range 



of sizes in the cellular context. 
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